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1 Introduction 



There are many powerful and systematically developed techniques in asymptotic theory for orthogonal 
polynomials. For instance, the steepest-descent method for integrals [9], the WKB (Liouville-Green) ap- 
proximation for differential equations [7], the Deift-Zhou's nonlinear steepest-descent method for Riemann- 
Hilbert problems [21 [3], and etc. Here, we intend to develop a simple, and yet systematic approach to 
derive asymptotic formulas for orthogonal polynomials by using their recurrence relations. Let {7r„,(x)}^o 
be a system of monic polynomials satisfying the recurrence relation 

■Kn+l{x) = {X - an)TTn{n) - bnlTn-lix), U > 1, (1.1) 

and the initial conditions 7ro(x) = 1 and vri(x) = x — oq. Note that for the sake of convenience, we have 
normalized the polynomials to be monic. To construct the asymptotic formulas of 7r^(x), we first set 

n 

■n:n{x) = Y\_Wk{x). (1.2) 

k=l 

It is readily seen from (II. ip that wi{x) = x — ao and 

Wk+i{x) = x - ak - —, k>l. (1.3) 

Wk 

When x is away from the oscillatory region of the orthogonal polynomials, it is easy to find an asymptotic 
formula for Wk{x) from (jl.Sp . Then, as we shall see, the asymptotic behavior of vr„(x) for x away from 
the oscillatory region can be obtained readily. When x is near the oscillatory region, we use a method 
similar to that given in [8l [10] to derive asymptotic formulas for general solutions of (jl.ip . The asymptotic 
formula of TTn{x) for x near the oscillatory region is then obtained by doing a matching. In the subsequent 
three sections, we will consider the following three cases: 

Case 1: = and bn = v? /{A-n? — 1). This case is related to the Legendre polynomials. 
Case 2: = and 6„ = n/2. This case is related to the Hermite polynomials. 

Case 3: = n"^ and = 1/4. This case was recently brought to our attention by M. E. H. Ismail. 

For simplicity, we use the same notations in the following three sections. Because each section is inde- 
pendent and self-contained, this will not lead to any confusion. 

2 Case 1: the Legendre polynomials 

The Legendre polynomials can be defined as [6l (1.8.57)] 

They satisfy the recurrence relation (6l (1.8.59)] 

(2n + l)xPn{x) = {n + l)P„+i(x) + ?iP„_i(x). 

For convenience, we normalize the Legendre polynomials to be monic. Put 

2"n' 

T^n{x) := 7 — — — P„(x). 

[n + l)n 

The monic Legendre polynomials {7r„(x)}^Q satisfy [6l (1.8.60)] 

7rn+l(x) = XTTnix) - _ _^ 7rn~l(x), n > 1, (2.1) 

7ro(x) = l, 7ri(x) = x. (2.2) 
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Theorem 2.1. As n ^ oo, we have 



for X in the complex plane bounded away from [—1, 1]. 
Proof. Set 



= n^'^(^)- (2-4) 



A;=l 

Prom (O), (12:2]) and (ITil). it follows that 



^,+i(x) = x-^A__L_, fc>i, (2.5) 
= X. (2-6) 

As /c — )• 00, we have 

, , X + a/x^ - 1 
w^fcW ~ ^ 

for X G C \ [—1, 1]. Here the square root takes its principle value so that Va;^ — 1 ~ x as x — )■ 00. Define 

»W := --±^^ (2.7) 



and 



:= (2.8) 



It is easily seen from ((23]), dM]) and (HSj) that 



ni(x) = ^. (2.10) 
ui x) 



We make a change of variable 



t = t{x) := (x - Vx^ - if. (2.11) 

It follows from (f277l) and (f2lT]l that 

u;(x)2 = l ^ = l + t. (2.12) 

Hence, the equations (|2.9p and (12.10p can be written as 

Ak'^t 1 

Wx) = l + t-^^^, fc>l, (2.13) 

u^(x) = l + i. (2.14) 
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Define Qo{t) := 1 and 

n 

Qn{t):=Y[uk{x), n>l. (2.15) 
k=l 

Prom ([233]) . (f2lH) and (f2J5]) . we obtain Q^{t) = l + t and 

Qn+l{t) = (1 + t)Q„,(i) - _^ Qn-l{t). 

From tliis recurrence relation, one can construct a generating function from which it is easily deducible 
that Qn{t) has the explicit expression 

fr-„ il(«-j + i/2), 

A simpler verification of this identity is by induction. Using the Lebesgue dominated convergence theorem, 
it can be readily shown that 

Q„(t)^(l-t)-i/2 

as n DO. Note that by ([23]), I^Bj and (12^51) . 7r„(x) = w{x)'^Qn{t). Thus, it foUows that 

TTnix) W{xni - t)-'/^ 

as n oo. This together with (|2?7l) and (|2TT]) yields ([23]). ■ 

Theorem 2.2. Lei 5 > &e an?/ /ixed small number. For x in a small complex neighborhood of the 
interval [—1 + 6,1 — 5], we have 



T^nix) ~ — 



l + sin6'\^/^ . /l-sin«^^/^ 



cos nO : — - — + sm n 



smy / V sm 



(2.16) 



as n ^ DO, where 9 = 0{x) := arccosx with < HeO < vr. 

Proof. To put the difference equation (j2.1|) in the form suggested by Wang and Wong [HI (2.1)], we let 

2"r(n/2 + l/4)r(n/2 + 3/4) 
Pnixj.- [r(n/2 + l/2)]2 ^^-^^^ 



From (I2.17p . it is easily seen that 

[r(n/2 + l/2)]2(n/2 + l/4) 
[r(n/2 + 1)]2 



2xpn(x) =Pn+l{x) +Pn-l{x). (2.18) 



Motivated by the form of the normal (series) solutions to second-order difference equations (see [101 (1.5)]), 
we assume 

Pn{x) ~ n°^\r{x)]^ {f {x) cos[n(p{xy\ + g{x) sin[n(^(x)] (2-19) 

as n — 7- oo, where r[x) and '^{x) are real-valued functions, whereas /(x) and g{x) can be complex-valued. 
We now proceed to determine the constant a and the functions r(x), V9(x), f{x) and g{x) in (|2.19|) . It 
can be easily shown from (j2.19p that 

Pn±i{x) ~ n"r"^"'^[(/ cos 93 ± (7 sin Lp) cos{nLp) + [g cos f + f sin (^3) sin(n(/9)]. (2.20) 
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Furthermore, by the asymptotic formula for the ratio of Gamma functions [U (6.1.47)] we have 

[r(n/2 + l/2)]2(n/2 + l/4) .^^.-2. 

[r(n/2 + 1)]^ = l + 0(n ) (2.21) 

as n ^ 00. Applying ^7l9\\ . and I^TTlh to ([2J8]) gives 

2x[f cos{nip) + g sin(n99)] ~ r[(/ cos (f + g sin 99) cos(n93) + (g cos if — f sin 99) sin(n(^)] 

+r^^[{f cos If — g sin if) cos{nip) + [gcosip + / sin 93) sm{n(p)]. 

Comparing the coefficients of cos(?7-99) and sm{n(p) on both sides of the last formula yields 

2xf = (r + r~^)f cos (f + (r — r^^)g sm(f; 
2xg = — (r — r~^)/sin93 + (r + r~^)g cos ip. 

Thus, we obtain from the above two equations 

X = cosh(log r) cos tp, = sinh(log r) sin (p. 

It can be easily seen that the only solution to these equations is logr = and p = arccosx. Recall that 

= e{x) := arccos x, < Re 6^ < tt. (2.22) 

Hence, we conclude that 

r = 1 and p = e. (2.23) 
Next we are going to determine the constant a in (j2.19p . Applying (j2.23p to (j2.19p gives 

Pnix) ~ n''[/cos(n6l) + gsm{ne)], (2.24) 

and 

Pn±iix) ~ f 1 ± -) [(/ COS 6* ± 5 sin 61) cos(n6l) + (g COS 6* =F / sin 9) sin(n6')] . (2.25) 



A combination of (f2l8]) . ([221]), and K2E\\ yields 

a 

22;[/cos(n0) + gsiii(n9)] ~ 2f cosO cos(n6) + 25fcos0sin(n0) H — [2(7 sin cos(n0) — 2/sin^sin(n0)]. 

n 

In view of (|2.22p . we obtain by matching the coefficients in the last formula 

ag sin 9 = 0, af sin 9 = 0. 

These equations hold for all x in a small complex neighborhood of [— 1 + (5, 1 — S]. Since / and g can not 
be identically zero, it follows that 

a = 0. (2.26) 

Thus, we have from (I2l9]l . (12:23]) and (12:26]) 

Pn{x) ^ f cos n9 + gsuin9 (2.27) 

as n — 7- 00. This formula holds uniformly for x in a small complex neighborhood of [—1 + 5,1 — 5\. 
Moreover, it follows from p.3p and (I2.17P that 

(r^ — \ 



for complex x bounded away from [—1,1]. Our last step is to determine the coefficients / and g in (I2.27P 
by matching the above two formulas in an overlapping region. With 6 and x given by (j2.22p . it can be 
shown that for Imx > 0, we have Im^ < 0. Thus, (]2.27p implies 



Pn{x) 



Z _L A ) pine 
2 ^ 2i ' 



Meanwhile, in view of x = cos e and \/x2 - 1 = isin6l by ([2:22]) . we obtain from ([2:28]) that 



Pn{x) ~ e 



ind 



Coupling the last two formulas gives 



,i(e-7r/2)"l^/^ 

2 sin 6* 

gi(9/2-7r/4) 



2 2i (2sin^)i/2- 
Similarly, matching ()2.27p with (|2.28p in the region Imx < yields 



/ 9 



-'i(6»/2-7r/4) 



2 2i (2sine)V2 ' 
From the last two equations of / and g we have 



/ 



1 + sini 



sin( 



1/2 



1 — sin ( 



sm ( 



1/2 



This together with ([2T7l) and i!a?Tl\ implies (1236]) . ■ 

3 Case 2: the Hermite polynomials 

The Hermite polynomials can be defined as [G] (1.13.1)] 

H.(.) = (2.r.F,(-"/2'-("-l)/2|-J5 

They satisfy the recurrence relation [G] (1.13.3)] 

2xHn{x) = Hn+iix) + 2nHn~i{x). 
For convenience, we normalize the Hermite polynomials to be monic, and put 

Tln{x) := 2-^Hn{x). 

The monic Hermite polynomials {'^n{^)}^=o satisfy [61 (1.13.4)] 

7r„+i(x) = XTTnix) - -7r„_i(x), u > 1, 
7ro(x) = 1, 7ri(x)=x. 

Theorem 3.1. As n ^ oo, we have 



7rniV2ny) ~ ( — j exp 



{n[y^ - yVy^ - 1 + iog(2/ + Vy^ - i)]} 



1/2 



(3.1) 
(3.2) 



(3.3) 



for complex y bounded away from the interval [—1,1] 
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Proof. Set 



(3.4) 



k=l 



It follows from ()3.1|) and (I3.2p that wi{x) = x and 



2wk{x) 



Let X — Xn ■ — yiny with y G C \ [— 1, 1]. It can be proved by induction that for real y and y ^ [—1,1], 
we have 



+ ^xl - 2k 



1 + 



1 ^Xfi 



2(4 - 2k) 8(4 - 2A:)5/2 
for all A; = 1, • • • , n. From these inequalities, it follows that 

Xfi + y4--2A; 



< Wk{Xn) < 



+ v^x2 - 2k 



1 + 



1 



2(4 -2fc) 



Wk{Xn) 



1 + 



1 



+ 0(n- 



(3.5) 



2{xl - 2k) 

as n — )• oo, uniformly in A; = 1, • • • , n. By using a continuity argument, it can be shown that the validity 
of this asymptotic formula can be extended to complex y G C \ [—1, 1]. Recall that x„ = y/2ny. By the 
trapezoidal rule 

/(i)-/(o) 



- E fik/n) = C f{t)dt + M_M + o{n 
n ^ Jq 2n 



we have 



l^^^^ x. + 74^ = -Elog(|/ + y^^^)-log(y + v/^) 

1 , y+^/y^-l 



n 



k=l 



k=l 

1 



log(y + Vy"^ - t)dt + — log 



2y 



log(y + y/y^ -1) 



/ -1/2-yvV^+^log 



(3.6) 



and 



k=l 



1 + 



^ 1 _ 1^ 1 

^2(x2-2/c) "n^4(2/2-fc/ 



2(4 - 2/c) 

/•^ dt ^ 1 

as n — )• oo. Applying the last two formulas and (|3.5p to (j3.4p yields 



n 



(3.7) 



n + x2 - 2n 



fc=l 

n\n/2 



n 

fc=l 



+ V4 - 2fc 
+ Y^4 ~ 2n 



n 

k=l 



1 + 



1 



y + V2?-l I exp n(y^ - 1/2 - y^/y^^^) 



2(4 -2/c). 



2y 



2/2-1 



n \"/2 



~ ) exp <i n 



[y2 _ y^fy^X + log(y + 7^)]} , 



thus proving (|3.3p . 
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Theorem 3.2. Let 6 > be any fixed small number. For y in a small complex neighborhood of [—1 + 
6,1 — 6], we have 



^2eJ (l-y2)i/4 
as n ^ oo, where 9 = 9{y) := arccosy 



cos 



n{6 — sin 9 cos 9) + — 



+ sin 



n{9 — sm 9 cos 9) + — 



(3.8) 



To prove the above theorem, we need a lemma analogous to Lemma 1 in [8] . For convenience, we use 
the notation 



y± 



\ 1/2 Q 

n \ y 

nil/ ^ ^ ^ 2n ^ 8n2 ■ 



(3.9) 



Lemma 3.3. Let (p{y) be any analytic function in a small complex neighborhood of [6, 1—6], we have 

cos[{n ziz l)ip{y±)] ~ cos(ny3) ^cos A =F — sin =F sin(n(^) ^sin A lb — cos (3.10) 

and 

sm[{n ziz l)ip{y±)] ~ sin(n(/?) ( cos A =F — sin A j ± cos(n93) ( sin A it — cos A j (3-11) 
as n ^ oo, where 



A = A(y) := ip{y) 



yf'iy) 



and 



^ = Ky) 



y^'jy) ^ y^v"{y) 



(3.12) 



(3.13) 



Proof. From (j3.9l) we have 

{n±l)ip{y±) ~ n{l±-)ip{y^ — + —) 

- n{ip±- + 
n 

where if denotes ^p{y), and A and /i are given in (j3.12p and (j3.13p . It then follows that 

cos{nip) cos(A lb fj./n) =F sin(n(^) sin(A it fj./n) 



cos[(n± l)ip{y±)] 
sin[(n± l)ip{y±)] 

This proves the lemma. ■ 
Proof of Theorem \3.SX Define 



~ cos(n(/?) ^cos A =F — sin A^ =F sin(n(^) ^sin A ± — cos A ) ; 
~ sin(n(^) cos(A it fi/n) it cos(n(^) sin(A it fi/n) 
~ sin(n(^) ^cos A =F — sin A^ it cos(n(^) ^sin A ± — cos A ) . 



Pn{x) := [r(n/2 + l/2)]-V„(x) 



(3.14) 
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We make a change of variable x = '■= \/2ny. It is easily seen from (|3.ip and (I3.14p that 
r(n/2 + l/2)\/2^ 



/o , i\ ypn{V2ny) = pn+i{V2ny) +pn-i{V2ny). (3.15) 
1 (n/2 + 1) 

As in (|2.19p . we now assume 

PniV2ny) ~ n°'[r{y)r{f{y) cos[n(f{y)] + g{y) sin[nv9(y)]} (3.16) 

as n — 7- oo. First, we shall determine the constant a and the functions r{y) and ^{y) in ()3.16|) . From 
(13:9]) and (l3lB . we have 

(V2(n±l)y±) 

~ (n ± l)"[r(y±)]"±H/(y±) cos[(n ± l)<^(2/±)] + ^(yi) sin[(n ± l)v^(y±)]}. (3.17) 



Moreover, it can be shown from (13. 9p that 

where r = r{y). Applying ([33]), (f3li)|) . (f3lT]) and (IXTHIl to (IXTTl) yields 

Pn±i{V2ny) ~ n°r"^"^e^^^ '^^''[(Z cos A ± (/sin A) cos(n99) + {g cos A =F / sin A) sin(n(^)]. (3.19) 
Here / and g stand for f{y) and g{y). By Stirling's formula [!_, (6.1.37)] we have 



r(n/2 + 1/2)VV2 1 



r(n/2 + 1) 4n 
A combination of (l3l^ . (I3l^ . ([3l^ and ([3:20]) implies 

2y[/ cos{n(p) + g sm{nip)] ~ re~^^ [(/ cos A + (7 sin A) cos{nip) + (51 cos A — / sin A) sm{nip)] 

+r~^e^'^ /^''[(/cos A — (/sin A) cos(n(^) + {gcosX + / sin A) sin(n(/?)]. 

Comparing the coefficients of cos(?i(/3) and sm.{n'^) on both sides of the last formula gives 

2yf = re-3'^'/2'^(/ cos A + 5 sin A) + r-^e3'^'/2'^(/ cos A - 5 sin A); 
2yg = re-y'''''^''{g cos A - / sin A) + e^"' Z"^'' {g cos A + / sin A). 

Solving the above two equations, we obtain 



yr'\ ■ 1 A y^' 



cosh I log r cos X = y, sinh ( log r sin A = 0. 

\ 2r J \ 2r J 

A solution is 

yr' 

logr = 0, cosA = y. (3.21) 

2r 

The first equation in (|3.2ip implies 

r = e'y^ (3.22) 
for some constant c S C. From (|3.12p and (|3.2ip . we have 

(f = ib(arccos y — y\/ 1 — y^) + dy^ + 2k'n 



for some constant c' G C and A; G N. Without loss of generality, we may take c' = and k = 0. Hence, 

(p = arccos y — y\/l — y^. (3.23) 
Next, we are going to determine the functions / and g in (I3.16p . From (13. 9p and (I3.22|) we obtain 

[r{y±)r^' = HyT (3.24) 

and 

/(y±)~/T^, giy±)^gT^, (3.25) 

where / = f{y) and g = g{y). Applying (f3l0]) . (f3lT]) . <!^?2^ and (l3:25D to (I3l71) yields 
Pn±i(\/2ny) 



(/ COS A lb (j( sin A) cos{nip) + ((j( cos A =F / sin A) sin(n(/9) 

cos(nv3) / .^f..^yf' , , . ^ , . yg' ■ . 

±aj cos fxj sm A =F -— cos A + at^ sm X + cos A — sm A 



n V 2 "^"^ 2 

, ^™(^'^) / 1 \ ^ ■ \^y9' X ^ • \ / \ , ■ \\ 

H I ±a(7 cos A =F l-i-g sm A =F -ir- cos A — aj sm A — /.i j cos A H — — sm A I . 

(3.26) 

A combination of (l3l^ . ([3l^ . ([CT]) and (f3:26]l gives 

1 — — ) [2y/ cos(n(/?) + 2y(7 sin(n99)] ~ 2/ cos A cos(n99) + 2g cos A sin(n(^) 



4n 



C0S(?1V9) ^ / X N 

- \lag sm A + 2;U(7 cos X — yg sm A) 



+^^^^^^(-2a/sinA- 2^/ cos A + y/' sin A). 

n 

In view of the second equation in (j3.2ip , we obtain by matching the coefficients in the last formula 

f 

— cos A + 2ag sin A + 2^g cos A — yg sin A = 0; (3.27) 
^ cos A — 2a/ sin A — 2^/ cos A + yf' sin A = 0. (3.28) 

Note from (f3l^ . (f3l3]) and (13:23]) that 

7/ 

A = arccos y, fi 



4^1 - y^ 

Hence, equations (|3.27l) and ()3.28p can be written as 



AagyT-y^ 

y ' Vi-y' 



/ + ^"^^^"^ + ^IL^ - = 0; (3.29) 



y 



(3.30) 



Set 

10 



We then have from ([329]), ([CT]) and (l3:3T]) 



u = ; v' = — 7^=. (3.32) 

Define 

6 = ^(y) := arccosy. (3.33) 
The solution of the system (|3.32|) is given by 

6 6 9 

u = Cicos- +C2sin-; t; = -Ci sin - + C2 cos -, (3.34) 

where Ci G C and C2 G C are two arbitrary constants. Consequently, we obtain from (j3.3ip that 

f = (i_^^2)i/4 (gi COS ^ + Casing ; (3.35) 

y2a 



(1 _ y2)i/4 ^ ^1 ^i"^ 2 + ^2 2 J ■ ^^-^^^ 

Applying (l3:22D . (I335D and (I336D to (f3l6]) yields 

PniV2^y) ~ n"e""J''y2"(l - y2)~i/^[Ci cos(n(/5 + 0/2) + C2 sin(nv? + 0/2)]. (3.37) 

This formula holds uniformly for y in a small complex neighborhood of [—1 + 5, 1 — 5]. Moreover, it follows 
from ([331) and (IXTil) that 

-\ 1/2 



p„(^/2;^y)~ -^exp{n[y2-y7^^ + log(y + 7^^)]} P +^^y_M (3.38) 



for complex y bounded away from [—1, 1]. Finally, we match the above two formulas in an overlapping 
region to determine the constants a, c, Ci and C2 in (j3.37p . For Imy > 0, it follows from ()3.33p 
that Im0 < 0; see a similar statement following (j2.28p . Furthermore, it can be shown from (I3.23P that if 
Imy > 0, then we also have Im^j < 0. [To do this, one first notes that ip'{y) is negative for y E [—1+5, 1—5]. 
Then, by the continuity of tp' , one concludes that Re (p'{y) < for y in a neighborhood of [— 1 + 5, 1 — 5] in 
the complex plane. Finally, the mean value theorem ensures that there exists a real number ^ G (0, Im y) 
such that ip{y) = ip(Rey) + i{lmy)ip' (Rey + i^), from which one obtains lniip{y) < 0.] Thus, ()3.37p 
implies 

p„(V^y) ~ n"e-^^^"(l - y'y^^ + e^^+'^Z^ 

Meanwhile, we have from (j3.33p and (j3.38p 

^ r ,7fl_^/o^ 1 ^/'^ 

Pn{V2ny) ~ — exp s n[y^ — iy\/ 1 — y^ + i arccos y] > 



Thus, we obtain from (j3.23p and the above two formulas that a = 0, c = 1 and 

Ci C2_ e-'^l^ 

T 27 ~ 2^ ■ 

Similarly, matching ()3.37p with ()3.38p in the region Imy < yields again a = 0, c = 1 and the equation 

CTi _ _ 

T ~ 2i" ~ V^' 
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Coupling the last two equations gives 

Ci = C2 



v/27r 

Therefore, we conclude that 

a = 0, c = 1, Ci = C2 ^ 



/27r 

This together with (|3l^ . (l3:23]l and (fSTSTl) yields 



4 Case 3: an open problem 

Recently, M. E. H. Ismail proposed the problem of finding asymptotic formulas for the orthogonal poly- 
nomials determined by 

7r„+i(x) = (x - n^)7r„(x) - ^7r.„_i(x), n > 1, (4.1) 

7ro(x) = 1, 7ri(2;)=x; (4.2) 

see O §6] and [H p. 370]. We first present a result for x not in the interval of oscillation. 
Theorem 4.1. As n — )■ oo, we have 

vr„(n2y) ~ (^-j exp + 1) log(^ + 1) - ( " 1) log( " 1)]} ( ^ ) (4-3) 

for complex y bounded away from [0, 1] . 
Proof. Set 



n 

■Kn{x) = 

k=l 



Uwkix). (4.4) 

k-- 

It follows from ()4.1|) and ()4.2p that wi{x) = x and 



u;fc+i(x) = x - k^ 



Awk{x) 

Let x = Xn ■= n?y with y E C \ [0, 1]. As with the case of Hermite polynomials, it can be shown that for 
real x and x ^ [0,n^], we have 

X - {k-lf -1 < Wk{x) <x - {k-lf + 1 

for all A; = 1, • • • , n. Thus, 

2k 2 Wk(x) 2k 
1 H < < IH . 

Consequently, 

2k 



I 2 \ 2 1 

Wk{n y) = n \y ^ 



(4.5) 



as n — )• oo, uniformly in A; = 1, • • • , n. By using a continuity argument, it can be shown that the validity 
of this asymptotic formula can be extended to complex y G C \ [0, 1]. In view of the trapezoidal rule 



n 1 

-Y^f{k/n)^ / f{t)dt + 



/(i)-/(o) 

2n ' 
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we have 

U2X ,1 1, y-1 



sr^ / k \ 1 
2^ log i^y - ^ J ^ ^ l°g(y - ^^)^* + 2 log ■ 



[(Vy + 1) iog(Vy + 1) - ( - 1) iog(Vy - i) - 2] + - log ■ 



y"io (^1+ '^^ \ V dt-io 



2 ° y 

and 

-dt = loff 

y-1 

as n — >• oo. Applying the last two formulas and (|4.5p to (j4.4p gives (|4.3p . ■ 

Next we give a result for 2; inside the interval of oscillation. 

Theorem 4.2. Let 6 > be any fixed small number. For y in a small neighborhood of [5, 1 — 5] in the 
complex plane, we have 



^nin'y) ~ (-1)^-12 sin(nvrVy) (^)'" {j^J^ y'^'i^ " y)""'^' (4-6) 



as n —)• 00 . 



To prove the above theorem, we will need a lemma analogous to Lemma 1 in [8]. As in p.9p . for 
convenience we set 

y±-=(^) y-yT — + %- (4.7) 

Lemma 4.3. Let ip{y) be any analytic function in a small neighborhood of [5, 1 — 5] in the complex plane, 
we have 

cos[(n lb l)(p{y±)] ~ cos(n99) ^cos A =F — sin =F sin(n93) ^sin A ± — cos A^ (4-8) 

and 

sin[(n lb l)(^(y-l-)] ~ sin(n(/?) ^cos A =F — sin A^ ib cos(n(/?) ^sin A ib — cos A^ (4-9) 
as n ^ CO, where 

\ = \{y):=^{y)-2y^'{y), (4.10) 

and 

11 = li{y):=y^'{y) + 2y\"{y). (4.11) 

Proof. From (j4.7p we have 

1 2y 3y 

{n±l)ip{y±) ~ n{l±-)ip{y=f b ^) 

1 2y^'^y^' 2y^^" 

~ n(lib- v'T + ^- + —) 

n n n'^ 

A /i 

n 
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where A and fx are given in (I4.10p and (14. lip . It then follows that 

cos[{n ± l)ip{y±)] ~ cos(n(/?) cos(A lb ^/n) =F sin(n(^) sin(A lb ^/n) 

~ cos(n(/?) (cos A =F ~ sin A ] =F sm{nip) (sin A ib — cos A 

\ Til / \ 

sin[(n lb ~ sin(n(^) cos(A lb /i/n) lb cos(n(^) sin(A lb ^/n) 

~ sin(n99) (^cos A =F — sin A ) ib cos(n99) (^sin A ib — cos A 

This proves the lemma. ■ 
Proof of Theorem \4-^ Define 



Pn{x) := \-^TTn{x). (4.12) 



We make a change of variable x = := n y. It is readily seen from (j4.ip and (j4.12p that 

(1 - y)Pn{n^y) = Pn+l{n^y) + ^^2^^_ ^-^2 Pn-li^^y)- (4-13) 

As in (|3.16p . we first assume 

Pn{n^y) ~ rf[r{y)T{f{y) cos[nv3(y)] + g{y) sm[nLp{y)\] (4.14) 

as n — 7- oo, and then determine the constant a and the functions r{y), f{y), g{y) and (p{y) in the formula. 
From (|47ll and (I4.14P we have 

Pn±i{n^y) = Pn±i{{n±lfy±) 

~ (n ± l)"[r(y±)]"±H/(y±) cos[(n ± l)v{y±)] + g{y±) sin[(n ± l)<^(y±)]}. (4.15) 



Moreover, it can be shown from (14. 7p that as n — )• oo, we also have 

where r stands for r(y). Applying (fO]) . (|OD . (1^91) and (Ii36]) to (|il3]l yields 

Pn±i{n^y) ~ n''r"=^^e'F2j/'''A[(jcos A lb 5 sin A) cos(nv3) + (^cos A=F / sin A) sin(n(/?)]. (4.17) 

A combination of (|4.13p . (|4.14p and (|4.17p gives 

(1 — y)[/ cos(n(^) + sin(n93)] ~ re"^^'^' /''[(/ cos A + gs'm A) cos(n(^) + ((7 cos A — / sin A) sin(n93)]. 

By comparing the coefficients of cos(n93) and sin(n(^) on both sides of the last formula, we obtain 

(l-y)/ = re-2^'>(/cosA + 5sinA); 
{l-y)g = re-23'^'/^(5cosA- / sinA). 

Thus, we have from the above equations 

{l-y) = re-^y'-'l'' cos A, = re-'^^'' I' sin A. 

The only solution is A = 0, and 

^g-2yr'/r = 1 _ (4.18) 
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With A = 0, we obtain from (|4.10p 



if = 



(4.19) 



for some constant c S C. Let R{y) := logr(y). From (j4.18|) . it is easily seen that R{y) satisfies a first-order 
linear inhomogeneous equation, whose solution is given by 



R{y) 



2^ 



log(l - s)ds 



Upon integration by parts, followed by a change of variable u = s^^^, one obtains 

R{y) = log(l -y) + 2yi/2j^rctanh^+ c'^ 
for some constant c' G C. Taking exponential on both sides gives 

r{y) = il-y) {l^y^'-^- 
Without loss of generality, we may assume c' = 0. Hence, 

.fe) = (i-,)(i±J)^. 

Next we determine the functions / and g in (|4.14p . From (|4.7p . (I4.18P and (I4.20p we have 

nil /I ..Mblr„/,.Mn i , U 



[r{y±)] 



1 + 



n(l - y) 

Furthermore, it is easily seen from ()4.7p and (j4.19p that {n ± l)ip{y±) ~ nip{y) and 



(4.20) 



(4.21) 



f{y±) ~ fiy) T 



2yf'iy) 



9{y±) ~ g{y) =f 



•^yg'iy) 



n n 
Applying the above formulas for functions r, ip, f and g to (I4.15P yields 



Pn±i{n'y)^n^r^{l-y)^Ul± 



a 



n 



1 + 



y 



n{l - y) 



2yf 



n 



cos{rup) + 5 =F 



2yg' 



n 



sm{rup) 



[One can also obtain this result from Lemma 14.31 since A = /.i = by ()4.19p .] This together with (14.131) 
and (I4.14P implies 



/ cos{nip) + g sin{nip) 



/ + - 



1 f yf 



+ af- 2yf 



cos{nip) + 



n \1 — y 

Comparing the coefficients on both sides of the last formula gives 



9+- 



( yg 



n \1 — y 



+ ag- 2yg' 



sm{nip). 



yf 
i-y 



+ af- 2yf' = 0, 



yg 

i-y 



+ ag- 2yg' = 0. 



Hence, 



/ = Ciy"/2(1 - y)-V2, g = C2y"/2(1 - y)-i/2^ (4.22) 
where Ci G C and C2 G C are two arbitrary constants. Applying (|4.19p . (|4.20p and (|4.22p to (|4.14p yields 



p 



.„(n2y)~n-y"/2(l-y)"-V2 



[Ci cos{nc^/y) + C2 sin(nc^)]. 



(4.23) 
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This formula holds uniformly for y in a small neighborhood of [6, 1 — 5] in the complex plane. Moreover, 
it follows from (gSD and I^^J^ that 

Pnin^y) ~ exp {n[( + 1) log( + 1) - (V^ " 1) log(^/y - 1)]} j (4.24) 

for complex y bounded away from [0,1]. At the final stage, we match the last two formulas in an 
overlapping region to determine the constants a, c, Ci and C2 in (j4.23p . In view of the equalities 
exp(±mcY^) = cos(ncy^) it isin(ncy^) and 

(1 - y)" [t~^) = '^''P ^""[^^ + + 1) - ( - 1) iog(i - Vy)]} > 

formula (j4.23p can be written as 

Pnin^y) ~ nV/'(i-y)"'/'exp{n[(Vy + i)iog(Vy + i)-(^/y-i)log(i- Vy)]} 

^ _ e--cvy +(^ + ^] e-c^y] , (4,25) 



2 2i J \2 2i 

Meanwhile, it follows from (|4.24p that for \my > 0, we have 

/ \ 1/2 

Til I y \ 

Pnin^y) ~ — exp {n[{^ + 1) log(^ + 1) - (y^ " 1) log(l - y^)] - imr^ - i7r/2} ( ) 

A comparison of the above two asymptotic formulas shows that a = 1 and c = vr or c = — vr. Without loss 
of generality, we take c = vr. Note that the function exp(mcy^) = exp(m7ry^) is exponentially small, and 
hence negligible in the region Imy > 0. By matching the last two formulas one more time, and ignoring 
the exponentially small term, we have 

Ci_C2_ e^^^/^ 
T ~ 27 ~ 27r 

With a = 1 and c = vr, we match (|4.24p with (I4.25|) in the region Imy < to obtain the other equation 

Ci C2_ 
T 2i" ~ 27r 

Upon solving the last two equations, we obtain Ci = and C2 = — 1/vr. Therefore, we conclude that 

a = 1, c = vr, Ci = 0, C2 = — l/vr. 
Combining this with and (103]) gives (fO]) . ■ 
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